Ballistic Spin Injection and Detection in Fe/Semiconductor/Fe Junctions 
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We present ab initio calculations of the spin-dependent electronic transport in Fe/GaAs/Fe and 
Fe/ZnSe/Fe (001) junctions simulating the situation of a spin-injection experiment. We follow a 
ballistic Landauer-Biittiker approach for the calculation of the spin-dependent dc conductance in 
the linear-responce regime, in the limit of zero temperature. We show that the bulk band structure 
of the leads and of the semiconductor, and even more the electronic structure of a clean and abrupt 
interface, are responsible for a current polarisation and a magnetoresistance ratio of almost the 
ideal 100%, if the transport is ballistic. In particular we study the significance of the transmission 
resonances caused by the presence of two interfaces. 

PACS numbers: 72.25.Hg, 72.25.Mk, 73. 23. Ad 



I. INTRODUCTION 

The controlled spin-dependent electronic transport 
through magnetic/nonmagnetic heterostructures is a cen- 
tral issue in the rising field of spin electronics .Bel In some 
cases, such as spin valves or giant magnetoresistance de- 
vices, the basic-science discoveries have lead to techno- 
logical applications within less than a decade. In other 
cases, however, as e.g. spin injection into the conduction 
band of semiconductors (SC), much remains yet to be un- 
derstood and achieved, experimentally and theoretically. 

The interest in spin injection from a ferromagnetic 
(FM) material into a semiconductor has been largely 
motivated by the proposed, but not yet achieved, spin 
field-effect transistor of Datta and Das .el There have been 
many tries, with increasiM^gc^rtjiMjleiMnstrate that 
such a device kJeasibleMaQBTOMBBM It has been 
already shownoQ that electrons in the conduction band 
of semiconductors can travel long distances without loos- 
ing memory of their spin. In parallel, many attempts to 
achieve spin-polarised currents have been made. ThcjiSft 
of magnetic semiconductors as leads of the junctionOu El 
would be a possibility, but they have the drawback of 
low Curie temperature, and thus would not be applicable 
at room-temperature. On the other hand, the attempts 
to use metallic ferromagnetic contacts were at first 
promising. Efforts to use InAs-based contactsEErEiJO due 
to their nice properties of an abrupt interface and an 
ohmic transition have resulted in very low current po- 
larisation, which might sometimes even be attributed 
to stray-field Hall or magnetoresistance effects.E-3 Sev- 
eral theoretical approaches based on the spin-diffusion 
or the Boltzmann equation h^Yf .shed light on the be- 
haviour of such systems Recently Schmidt et 
al.il revealed a basic obstacle for succesful spin injec- 
tion, namely the conductivity mismatch between the FM 
and the SC, resulting in too low current polarisation un- 
less the FM contact is almost 100% spin polarised. Their 
conclusion holds in the diffusive regime, when one can use 
a resistor model for the FM/SC/FM heterostpucture. To 
overcome this fundamental difficulty, Rashbalia and Fert 
and JaffreseJ have proposed that the FM and SC parts 



should be separated by a tunneling spin-polarising slab, 
the high resistance of which would balance the drawback 
of the conductivity mismatch of a direct contact. In par- 
allel, and independently from thesfi-suggestions, there has 
been the observation of GrundlerEll that a ballistic tran- 
sistor would allow for a higher current polarisation than 
a diffusive one, and that this should be realisable if a 
two-dimensional electron gas was used. Already, simple 
model calculationsJ^Eiluj based on a free-electron ap- 
proach of the electronic structure of the leads, have shown 
that ballistic transport can give spin injection efficiencies 
of a few percent. Much more is seen, though, when one 
takes into account the full band structure of the FM ma- 
terial and the self-consistent electronic structure .of the 
interface. Indeed, as first proposed by KirczenowJ^J one 
can have ideal spin filters if the FM Fermi surface of only 
the one spin direction, when projected to the plane of the 
interface, has no states in the part of the two-dimensional 
Brillouin zone where the conduction band starts, so that 
there is no propagation into the SC from this spin chan- 
nel. This "selection rule" unfortunately does not apply in 
certain interesting systems such as Fe/GaAs or Fe/ZnSc. 
Nevertheless, as shown by Wunnicke et al.f3 in these sys- 
tems the interface reflectance is so much different for the 
two spin directions that one gets spin injection ratios as 
high as 99% in an ab initio ballistic calculation. Apart 
from the theoretical efforts, there arepSome very encour- 
aging recent experiments of Zhu et a/o giving already a 
few percent of current polarisation in Fe/GaAs(001). 

In the current article we present ab initio calculations 
of ballistic spin-dependent transport in Fe/GaAs/Fe and 
Fe/ZnSe/Fe trilayer heterostructures grown epitaxially 
in the <001> direction emulating a spin- valve geome= 
try. In this way we extend the work of Wunnicke et al.ti 
to include spin injection and detection. We show that 
the presence of the two spin-filtering interfaces increases 
the current polarisation even closer to the ideal 100%, 
and we also calculate the high magnetoresistance ratios 
of these structures, which is also approaching the ideal 
100%. We observe interesting interference effects due to 
the presence of two interfaces, and give an aspect of the 
whole problem that brings it in close connection with the 
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theory of ma gne tic tunnel junctions as it is described in 
Refs. ^5] and ]2(| Our results thus stress that epitaxial 
junctions operating as close as possible to the ballistic 
regime can form almost ideal spin filters and can exhibit 
extremely high magnetoresistance ratios. 

The article is organised as follows: in Section II we 
present the basic formulae of our ab initio approach. In 
Section III we describe the junctions to be calculated 
and the approximations made. The role of the symme- 
try of the wavefunctions in transport through Fc/SC/Fc 
junctions is explained in Section IV. Sections [v| and VI 
contain the results for the current spin polarisation and 
a discussion o f int eresting interference resonance effects, 
while Section VII is devoted on the case of antiparallcl 
orientation of the leads and the magnetoresistance prop- 
erties. Finally we discuss the limitations of our approach 
and conclude with a summary in Section VIII. 



II. METHOD OF CALCULATION 

Our calculations are based on density-functional the- 
ory in the local spin density approximation (LDA). We 
employ the screened Koxringa-Kohn-Rostoker (KKR) 
Green's function methodEH to calculate the electronic 
structure of the systems. In this multiple-scattering ap- 
proach, the one-electron retarded Green's function at en- 
ergy E is written in terms of local wavefunctions i?£ (r) 
and H2 (r) (regular and irregular solutions of the single- 
site Schrodinger equation, respectively, characterised by 
the angular momentum index L = (Z,m)), centered at 
lattice sites R„ and R n ', as 

G(R„ + r, R n . + r') 
= -n/£^i*2(r < ) J ff2'(r > )<W 



(1) 
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with G™l'(E) the so-called structural Green's function 
describing the intersite propagation; r< and r> are re- 
spectively the shorter and longer of r and r', and atomic 
units have been used (e = —y/2, % = 1, m = 1/2). 
The structural Green's function is related in turn to the 
known Green's function of a reference system via an alge- 
braic Dyson equation. For more details on this we refer 
the reader to Refs. ^8] and |2jj. 

The systems consist of two half-infinite (Fe) leads, as- 
sumed to have perfect periodicity otherwise. Sandwiched 
between these leads is an "interaction" region where a dif- 
ferent material (SC) can be placed and where the scat- 
tering of the Bloch waves takes place. The interaction 
region and the two leads have common in-plane Bravais 
vectors, i.e. in-plane (x-y) periodicity (perpendicular to 
the growth direction). If needed, larger (non-primitive) 



two-dimensional unit cells are taken to match the lattice 
constants of the materials. The two-dimensional period- 
icity of the layered systems allows to Fourier-transform 
the Green's function in the x and y directions, obtaining 
a two-dimensional Bloch vector k|| = (k x , k y ) as a good 
quantum number, and retaining an index i to charac- 
terise the layer in the direction of growth z. The Green's 
function connecting the layers i in the left lead and i' in 
the right lead is then written 

G(Ri + Xu + r, Ri' + Xw + r') 



47r 2 5sBZ Jsbz 
xJ2Ri(r)Gf L ,(k r ,E)RUr') 



(2) 
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where \v an d Xv' are in-plane lattice vectors, Ri is the 
interlayer lattice vector, SBZ is the surface Brillouin zone 
of the system and <Ssbz its area. In this equation each 
layer i is assumed to have a unique atom type, hence only 
the index i suffices to characterise the local wavefunction. 
In the case of more inequevalent atoms per layer, an ex- 
tra index is introduced to account for the propagation 
between different kinds of atoms. Moreover, in the case 
of ferromagnetism, the Green's function is different for 
each spin direction a =f or [. 

For the calculation of the conductance in lin- 
ear response all the information needed is contained 
in the Giles's function. In the Landauer-Biittiker 
approach, EllEH which identifies the ballistic conductance g 
with the transmission probability of the conducting chan- 
nels, one has 



■9= ^EE T ( k ll'^» 

<T,kii /X,// 



(3) 



relating the transmission probability T per channel to 
the conductance g. Here each channel is characterised 
by the band index /x, the k|| -vector and the spin a of the 
incoming electrons, and similarly by the primed indices 
for the outgoing electrons, both having the same Fermi 
energy Ep. Conservation of spin due to assumed absence 
of spin-orbit scattering, and of k| due to two-dimensional 
periodicity, have allowed us to omit the summation over 
a' and k'| in the outgoing electron channels. We follow 
here the formalism of Baranger and Stone,El relating g to 
the spatial derivative of the Green's function connecting 
a cross-sectional plane in the left lead (L) to one in the 
right lead (R). It is assumed that these planes lie in 
the asymptotic regime, where interface perturbations and 
evanescent interface states are no longer present. The 
formula for the k||-projected conductance g{k.\\,cr) per 
two-dimensional unit cell surface area and spin a reads 
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where the symbol d z stands for 

f{v)T z g{v) = f(v)d z g{v) - (d z f(r))g(r). (5) 

The conductance is evaluated only at the Fermi level Ep 
since we are at the limit of zero temperature. The com- 
plex conjugation in the last term of eq. (^) comes from 
conversion of the advanced Green's function to the re- 
tarded one by conjugation and exchange of r and r'. 
G(r, r'; kii) is given by the two-dimensional Fourier trans- 
form of (0) . By virtue of the Fourier transformation, the 
integration in (Q) is not performed over the whole lead 
cross-sectional area, but only over a two-dimensional unit 
cell. The total conductance per two-dimensional unit cell 
surface area for each spin channel is then 

g a = [ d 2 k l]g (k h a). (6) 

JSBZ 

Current conservation guarantees that the result is inde- 
pendent of the position of the cross-sectional planes of 
integration, as long as they are chosen in the asymptotic 
region. Details about the evaluation of the conductance 
will be given elsewhere£3 

The formula we use for the conductance has be pcoven 
to be equivalent to the Landauer-Biittiker formula.123 The 
conductance we calculate is then fully ballistic; diffuse 
scattering is assumed to be absent. Our approach also 
ignores spin-orbit scattering and any spin- flip events. We 
must also note that the semiconductor band gaps are 
known to be underestimated in the LDA by a factor of 
about 50%. This can have some quantitative significance, 
but the trends of our results are expected remain unal- 
tered even if we choose to enlarge the gap artificially. 

In the calculations, the atomic sphere approximation 
(ASA) for the potentials is used, i.e. they are assumed to 
be spherically symmetric around each atomic site and 
to occupy an atomic volume; on the other hand, the 
full charge density, rather than its spherically symmet- 
ric part, is taken into account. Moreover, we treat the 
systems nonrelativistically. An angular momentum cut- 
off of Z max = 2 has been taken for the wavefunctions and 
Green's funcions in the self-consistency procedure. 

III. THE SYSTEMS UNDER STUDY 

We study the spin-dependent transport through 
Fe/GaAs/Fe and Fe/ZnSe/Fe junctions. The junctions 
are supposed to have grown epitaxially on (001) inter- 
faces, and in an ideal way so that the transition from one 
material to another is abrupt. Absence of interdiffusion 
and of disorder are assumed; in this way, we are dealing 



,r^k||;£ F )^^G;(r,r';k||;£; F ) (4) 

I 

with a system grown in the z direction and being trans- 
lationally invariant in the x and y directions . The Fe 
leads are supposed to be infinite, while the semiconduc- 
tor thickness is varied from 41 to 97 monolayers (ML). 
In such thicknesses, the evanescent interface states in the 
semiconductor are expected to have decayed to insignif- 
icance compared with the Bloch wavefunctions, so the 
transport will be mediated through propagating states. 

Throughout the system, the experimental Fe lattice 
constant of ap c = 2.87lA is used. Thus, all atoms sit 
on ideal positions of an underlying bcc lattice. In par- 
ticular, in the SC part, the zincblende structure can be 
easily seen to fit on such a lattice, with half of the bcc 
sites occupied by Zn and Se (or Ga and As) atoms and 
the rest occupied by vacancies. Viewed in this way the 
consecutive positions of the atoms in the cubic diagonal 
of the bcc lattice are (Zn, Se, vacancy, vacancy). The 
zincblende lattice constant is then twice the one of the 
bcc. One can see, that using 2 x ap c = 5.742A in the 
SC part results only in a slight mismatch of less than 
2%, the experimental lattice constants being 5.654A for 
GaAs and 5.670A for ZnSe. In all cases, Zn termina- 
tion of the ZnSe spacer and Ga termination of the GaAs 
spacer was considered. As shown in Ref. |24|, the spin 
polarisation of the current through the single interface 
for the other terminations (Se and As) is also extremely 
high, and from the analysis of Sections [v|- |VIl| it follows 
that the two-interface junctions for those terminations 
will have qualitatively the same properties as the ones 
studied here. The two planes L and R used for the in- 
tegration were 6MLs away from the interfaces in the Fe 
region, where the asymptotic regime is assumed to have 
been reached. Variation of this distance causes insignifi- 
cant changes in the results. 

In a system as the ones we are considering, the Fermi 
level will be naturally determined by the infinitely long 
Fe leads. But in the spacer material, two or three mono- 
layers after the interface, the potentials and the charge 
density must be almost bulk-like. For this reason, the 
potentials of the inner atoms of the spacer will be auto- 
matically adjusted to the Fe Fermi level by a constant 
shift which is the result of the interface dipole layer. The 
self-consistent calculation of the potential close to the 
interface is then essential. 

Since we want to inject electrons into the SC conduc- 
tion band, we must emulate in some way a gate volt- 
age, or energy shift, acting on the SC potentials in order 
to lower the conduction band minimum slightly under 
the Fermi level. This artificial shift is different than the 
one just mentioned above, and it enters as a parameter 
in our calculations. We avoid disturbing the interface 
electronic structure, which is strongly influenced by the 
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metal-induced gap states, and proceed as follows. EJ The 
first two SC monolayers adjacent to the interface are kept 
as calculated by a self-consistent calculation of a 9ML- 
thick SC slab sandwiched between infinite Fe leads. The 
same applies also for the first neighbouring Fe MLs. Hav- 
ing saved the interface in this way, we take for the rest 
of the SC spacer (third up to last-but-two ML) the bulk- 
like potential that we find for the atoms in the middle 
of this Fe/9ML SC/Fe junction. This is justified, since 
it is known that the potential stabilises quickly as men- 
tioned previously. The emulation of the gate voltage is 
achieved by applying to this potential an extra shift such 
that the conduction band minimum E c of this bulk-like 
structure falls slightly under the Fermi level E F of the 
whole structure: 

E F = E C + E . (7) 

The parameter Eq, characterising the assumed gate 
voltage, is varied in our calculations over three val- 
ues: 20mRy, lOmRy, and 5mRy (272meV, 136meV and 
68meV, respectively). In this way we are able to view the 
approach to small values as a limiting procedure; as we 
shall see, these values are already in the limit of large spin 
polarisation of the current and magnetoresistance. View- 
ing the semiconductor part, the small values of Eq mean 
that the energy dispersion relation is nearly parabolic, 

m)~E c c-Lk^-^m + ki) (8) 

where m* = (d 2 E/dk 2 )~ 1 is the effective mass, and that 
the Fermi wavenumber k F is very small: 

E F -E c = E = -^k 2 F = ^ 7 (]J l +k 2 z ). (9) 
m m " 

These relations are relevant in the semiconductors consid- 
ered here because of their direct band gap at the T-point. 
Because Eq is very small, we have a very small Fermi 
sphere in the semiconductor. For this reason, very few 
channels kii will be able to conduct, namely those close 
to the center of the Brillouin zone with |k||| < k F . For 
the rest k z becomes imaginary and represents decaying 
wavefunctions. These can give rise to a tunneling cur- 
rent, but for the larger spacer thicknesses they are small 
compared to the contribution from the central part of the 
Brillouin zone. In any case they are always included in 
the calculation. 



IV. THE IMPORTANCE OF SYMMETRY 

As mentioned in the previous section, we are expecting 
contributions to the current only from the central part of 
the (001) Surface Brillouin Zone (SBZ), i.e. from k|| close 
to the T-point. In view of this, we will examine the ex- 
pected behaviour for states exactly at ky = 0, and argue, 
and show in fact in the calculations, that by continuity 
the close-by states will behave similarly. 



To begin with, we must clarify that the two- 
dimensional unit cell and the SBZ are determined by 
the SC part, since one SC lattice constant is assumed to 
match exactly two Fe lattice constants in our case. The 
states with k|| = can be examined in a great extent 
through their symmetry properties, since the z-axis re- 
mains invariant under many point-group operations. The 
single Fe (001) surface is characterised by the symmetry 
group Cztt,, having eight operations: a fourfold rotation 
axis (here the z-axis) plus reflections over the planes con- 
taining the z-axis and the xy diagonal or antidiagonal. 
But the zincblende (001) surface has the symmetry group 
C2V1 having four operations (a twofold rotation axis plus 
the reflections over the xy diagonal and antidiagonal), 
and being a subgroup of the former. As a result, the 
combined Fe/SC interface is characterised by the group 

The idea now, in view the Landauer approach and 
Eq. (||) , is to investigate the incoming states at the Fermi 
level deep in the Fe lead, as incoming channels, in order 
to see if their symmetry properties allow them to couple 
to SC propagating bulk states, and then see if these in 
turn are allowed (by symmetry) to couple to the outgoing 
states which propagate deep in the other Fe lead. The 
different character of the Fe states for majority and mi- 
nority spin will give us in this way hints about the spin 
polarisation of the current. This procedure can be used to 
propose theoretically ideal spin filter systems. But note 
that in this way we can only find which channels are ex- 
cluded from transmission by symmetry. As we shall see, 
some channels can be almost blocked for other reasons, 
contributing (by their absence from transmission) to the 
spin injection effect. 

We can now turn our attention to Fig. [l], where the 
energy bands of Fe, ZnSe and GaAs are drawn for k x = 
k y = in the /c z -direction, which is the one of inter- 
est as discussed earlier. Each of them is named by the 
irreducible representation to which it belongsEj for ro- 
tations around the A-axis (i.e. k z ). For example, the 
state labeled "1" corresponds to the Ai, which means 
that the states are invariant under all group operations 
(rotations around the z-axis); the label 2' refers to the 
A2' representation, being invariant under reflections from 
the planes containing the z-axis and either the xy diag- 
onal or antidiagonal. But we must note that for Fe the 
nomenclature refers to the C^ v group, while the sym- 
metry group of the whole system as well as of the bulk 
semiconductor is Ci v . Therefore we must use the com- 
patibility relations between the two groups, that show us 
which representations of C^ v have nonzero projection in 
each representation of Ci v . These can be found, for in- 
stance, in Ref. |34|. In our case we see that, at the Fermi 
level, only one band exists in the semiconductor (both 
for GaAs and ZnSe) , and it belongs to the representation 
Ai(C2u) (in parentheses we specify the point group to 
which the representation belongs). With this represen- 
tation, only the Ai(C4,j) and ^(C^) states of bulk Fe 
are compatible. This means that incident states of only 
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FIG. 1: Energy bands of bulk Fe (left) together with bulk 
ZnSe (centre) and bulk GaAs (right) along the A-direction 
(k z ), corresponding to F — H in bcc (Fe) and to T — X in fee 
(zincblende). For Fe, the black lines represent majority-spin 
states, and the gray lines minority-spin states. The potentials 
of GaAs and ZnSe have been appropriately shifted so that the 
Fermi level falls slightly in the conduction band. Each band 
is named by the corresponding irreducible representation of 
the point group; e.g. 1 means the Ai representation, 2' the 
A 2 < etc. For the notation see, e.g., Ref. B4L Note that the k z - 
axes at the semiconductor plots should actually be half the 
size shown, since the lattice constant is assumed double the 
one of Fe. Backfolded bands due to the doubling of the Fe 
two-dimensional unit cell are unimportant and not shown. 



these symmetries can couple to the semiconductor con- 
duction states (or even to each other, near the interface) 
and propagate into the SC spacer, while the rest, A2 {C^ v ) 
and A 5 (C4„), are totally reflected at the interface. 

Now, the energy bands of Fe can are quite different for 
majority vs. minority electrons near the Fermi level, due 
to the exchange splitting. At Ep the majority electrons 
have a Ai(C4t,)-state that can couple to the semiconduc- 
tor, while this is absent for the minority-spin carriers. 
For these, on the other hand, a A2' (C^-band exists 
that can do the job. We note in passing that this ab- 
sence of A 1 (C , 4 t ,) is due to the so-called s-d- hybridisation 
gap which splits this band in two and which happens to 
fall around Ep for the minority-spin states. 

If the A2' (C4t,)-band were absent, or if it could not cou- 
ple to the Ai(C2t,)-band of the semiconductor, we would 
be facing an ideal spin filter: only majority-spin would 
be able to propagate. Even in our case, however, we shall 
see that almost ideal spin filtering will occur, because the 
two kinds of states, Ai {C± v ) and A 2 ' (C*4„), have very dif- 
ferent transmission probabilities through the interfaces 
such that the A2<(C4t,) channels are nearly blocked. 
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FIG. 2: Spin-dependent conductance for Fe/ZnSe/Fe (top) 
and Fe/GaAs/Fe (bottom) junctions, as a function of k x , 
for the parallel magnetic configuration of the leads. The 
majority-spin conductance is illustrated in the left pannels 
and the minority in the right. Several SC spacer thicknesses 
are considered (49ML to 97ML), and gate energy shifts of 
Eq = 5, 10, and 20mRy. For the minority-spin case a magni- 
fication from 10 2 to 10 6 (see inset numbers) has been used to 
bring the graphs to the same scale. The fc^-values for ZnSe 
are 0.038, 0.056, and 0.083, and for GaAs 0.021, 0.031, and 
0.050 for Eq = 5, 10, and 20mRy, respectively. 



V. RESULTS FOR THE SPIN-DEPENDENT 
CONDUCTANCE 

The spin-dependent conductance as a function of k|| 
for several spacer thicknesses is shown in Fig. |^ for GaAs 
and ZnSe spacers, and for several energy shifts Eq. The 
wavevector kn has been taken along the T-X cubic direc- 
tion, which in the two-dimensional geometry corresponds 
to T-M. It is most convenient to express k|| in units of 
27r/asc = 27r/(2aF ), since this corresponds to the two- 
dimensional periodicity of the whole system; henceforth 
these units will be implied but omited for simplicity. The 
calculated values of kp are given in the caption of Fig. ||. 

The first evident observation is that the conductance 
practically vanishes for k» > hp, as expected. This effect 
shows up clearer for the thicker spacers. Thus we can 
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see that, as -Bo rises and the Fermi sphere in the SC 
becomes larger, the cutoff in conductance moves to higher 
values of fey, exactly as kp. As mentioned earlier, for 
larger values only evanescent states can exist, giving rise 
to a very small tunneling current which dies out as the 
spacer gets thicker. Nevertheless our calculations show 
that these states dominate the behaviour in the small 
thickness region. 

One can clearly see that the minority-spin conductance 
is lower by orders of magnitude than the majority coun- 
terpart. This is clearly the effect of the Fe minority A2'- 
state not being able to couple well with the SC Ai-state 
at the interface. The reason for this is that the A2' -state 
consists locally of d^-like site-centered orbitals. These 
point in-plane and are quite localised, so they cannot 
overlap very well with the SC Ai(C2„) orbitals. More- 
over, the SC Ai(C 2v ) band consists of s, p z , and d^-like 
states. The latter are in fact the ones that do couple 
to the A2' minority band of Fe. But we must note that 
such d^-like SC states are not inherent to the SC atoms, 
but rather induced as a distortion to the inherent sp SC 
orbitals by the neighbouring atoms sitting in the tctrahe- 
dral positions and giving a directional preference; in this 
sence they appear just as a correction when we use an an- 
gular momentum basis. As we depart from the T-point, 
other Fe minority orbitals (the continuations of the A5 
and A2 bands) begin to couple slowly, so the transmission 
increases. 

In contrast, the Ai(C4„)-band present in the majority- 
spin states consists locally of d z 2 , as well as s and p z - 
like atomic orbitals; these, pointing partly into the SC 
and being more extended, favour a better overlap and 
bonding with the SC states. Thus the reflectance of the 
interface is by far stronger for the minority-spin electrons, 
and a strongly polarised current results. 

The arguments presented here show that one needs a 
clean and abrupt interface, so that k|| is conserved. In 
the case of kn -violation due to diffuse scattering the ef- 
fect of spin selection will be reduced. Indeed, the to- 
tal (i.e. ky -integrated) density of states of Fe at Ep is 
higher for the minority-spin than for the majority-spin. 
On these grounds one would expect even a negative-cur- 
rent polarisation, in similarity with Julliere's modeS for 
spin-dependent tunneling; this might be the case if strong 
diffusive scattering intermixes the scattering ky -channels 
in a completely random way. Thus it is the specific selec- 
tion rule imposed by the interface in the ballistic regime 
that causes the strong positive current polarisation. We 
may also note that for other interfaces, such as (110) or 
(111), the symmetry of the various incident states is dif- 
ferent than in (001), and the selection rule might not be 
as strong; an ab initio calculation is necessary in order 
to judge this. 

Note that the same effect appears when one looks at 
tunneling, rather than spin injection, in these structures. 
This is demonstrated in Ref. where the tunneling 
through the semiconductor is also confined to the states 
close to T; there, the majority Fe state of symmetry Ai 



couples much better at the interface than the minority 
state of symmetry A2' , while both propagate with equal 
difficulty afterwards in the SC, as evident by the equal 
decay rate. In that case, of course, one must consider 
the complex band structure of ZnSc in the gap region 
as the analytical continuation of the conduction band pf 
the same symmetry for the interpretation of the effect,c3 
rather than the real conduction band structure, but both 
tunneling and spin injection can be viewed in this respect 
in a unified way. 

VI. INTERFACE REFLECTANCE AND 
QUANTUM WELL STATES 

Another interesting feature is the multi-peaked struc- 
ture of g{kn). This is an interference effect to the discus- 
sion of which we turn now. We start with the observation 
that the presence of two Fe/SC interfaces can give rise 
to interference effects due to the coherent multiple re- 
flection of the electrons between them. So, one expects 
resonances in the transmission, similarly to the case of 
a square barrier of finite length—met by free electrons of 
energy higher than the barrier £3 More concretely, let us 
assume that the transmission through each Fe/SC inter- 
face (1 or 2) has an amplitude ti t2 and the reflection r\_ 2 . 
These contain the phase shifts </>i.2, that the wavefunc- 
tion obtains for each reflection, plus a phase factor of 
e ik z D £ Qr wave propagation from side to side of the 
SC slab of thickness D leading to a phase of 2k z D for a 
come-and-go. A resonance in transmission will be formed 
whenever there is constructive interference after a num- 
ber of comes-and-goes of the wave; i.e. one has to sum 
up the series 

hot = hh + tir 2 rit 2 + tir 2 rir 2 rit 2 H 

= ll 1 - |r 1 ||r 2 |e i ( 2fc - £, +^ 1 +^)* 2 ^ 

in order to find the maxima in transmission. If the two in- 
terfaces are the same, as in Fe/SC/Fe with parallel mag- 
netic orientation of the leads, the single-interface proba- 
bilities T s j of transmission and i? s ; = 1 — T S j of reflection 
are equal for the two interfaces (the Fe/SC interface is 
equally hard to cross in either direction), and by squar- 
ing the previous equation one finds the total transmission 
probability to be 

Ttot = Ktotl 2 (11) 

T 2 

_ SI 

1 + (1 - T si ) 2 - 2(1 - To) cos(2k z D + fa + fa) 

where we have used the fact that T s \ = \t\t 2 \ is valid in 
this case; Eq. (jl^) is equivalent to the formula of Airy 
for a Fabry-Perot interferometer. This function is clearly 
oscillatory in k z D, and it exhibits a maximum of Xtot = 1, 
i.e. a resonance, whenever the condition for constructive 
interference is met: 

fa + fa + 2k z D = 2irn, (12) 
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with n an integer. 

For a given thickness D, variation of fen will cause 
variation of_fc z , and this will lead to these resonance 
phenomena.^ This is realised by combining eqs. JM) and 
(|l2|), so that the multi-peaked structure in Fig. gis ex- 
plained. To see what one expects qualitatively, we com- 
bine eqs. @ and (|l2|) to get 



Fe/ZnSe/Fe 



k F 



klD = 2wn 



(13) 



as a resonance condition. For zinc-blende structures, 
where k varies between and 1 in units of 2ir/asc, and 
o-sc = 4ML in the (001) direction, the condition relates 
feii to the number of monolayers A^ml: 



k\ - fcjj N M l = 2n + (0i + 2 )/vr. 



(14) 



Naturally, cf>i and 02 depend on k|| . This formula can be 
seen to give three resonances already for Nml = 100 and 
k F = 0.05(27r/a SC )- 

Between the maxima there are minima of T tot = 
Tj/(2 - T si ) 2 . For low values of T si the halfwidth of 
the resonance becomes very small; this is reflected at 
the minority-spin conductance where the resonances are 
much more narrow and peaked, with extremely low val- 
ued valleys between them, and thus their k|| -integrated 
contribution remains insignificant compared to the ma- 
jority one. These arguments also demonstrate that the 
interference effects arc in practice unable to invert the in- 
jected current polarisation, in contrast tp-.what has been 
predicted by recent model calculations .cfl We also note 
that, for E F -> E c , T si cx k z oc y/E F - E c s3 Then 
for a given spacer thickness T tot goes to zero linearly as 
T to t oc E F — Ec, while the first resonance appears for 
a thickness D ICS increasing to infinity as l/\/E F — Ec- 
In the model described here, one can readily substitute 
the values of T s \ from .a-single- interface calculation,! 2 .!! and 
get the correct trend£3 Nevertheless, in the calculations 
we cannot observe a perfect resonance of transmission 
one, because perfect coherence is destroyed by a very 
small but nonzero imaginary part of the energy, numeri- 
cally necessary for the calculation of the retarded Green's 
functionJ!3 

Another aspect of the matter is this: at the resonance 
values of k z D we have also a formation of quantum well- 
like states in the spacer. They are not bound, since the 
reflection is not total; the "interactive" change in the 
integrated density of states for each kii because of them, 
compared to the bulk Fe, is§3 

AN(E) = --Imln(l - |r 1 ||r 2 |e l(2fe ^ +<# ' 1+fe) ) (15) 

per spin direction. Whenever such a quantum well state 
is met, a resonance in the transmission probability is ex- 
pected; the larger | r*i | j r*2 j is, the more peaked and lo- 
calised in energy iSj-the change of the DOS and the trans- 
mission resonance.^ 
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FIG. 3: Majority- (left) and minority- (right) spin conduc- 
tance at k|| = as a function of the ZnSe spacer thickness. 
The oscillations of period 2n/(2kp) are evident; the values of 
2n/(2k F ) are 24.1ML, 35.7ML and 52.6ML for E = 20, 10 
and 5mRy, respectively. The peaks are much more violent 
for the minority-spin case (note the logarithmic scale there) 
because of the greater confinement due to stronger interface 
reflection. 



Dual to the oscillations of g in fe-space are oscillations 
in real space, when the spacer thickness D is varied while 
k|| is kept constant. As can be read out from eq. (|l2|), one 
expects a thickness period of 27r/ (2k z ), which for k|| = 
becomes 2ir/{2kF)- Indeed, in Figs. || and [I| we can see 
this oscillatory effect on the majority-spin conductance 
(left panels) for both ZnSe and GaAs spacers, with ex- 
actly the predicted period. The period gets longer for 
lower energy shifts, since they correspond to lower k F . 
On the other hand, larger fc|| will result in larger peri- 
ods, until the limit value of fe|| = k F ; after that k z be- 
comes imaginary, and one has attenuation rather than 
propagation of the wave, described by the complex band 
structure, as in a tunnel junction. 

Similarly, the minority-spin conductance oscillates 
with the same period as seen in Figs. |§| and ^ (right pan- 
nels) , but for the reasons mentioned before the peaks are 
much more pronounced; note that in this case a loga- 
rithmic scale was used for the intensities. It should be 
noted that there is, in particular for GaAs, an initial 
exponential decrease in the conductance, before the os- 
cillations start, as can be seen from the characteristic 
linear behaviour in the logarithmic scale. This origi- 
nates from decaying states with complex Bloch vectors, 
which contribute to the conductance by tunneling. In- 
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FIG. 4: Same as in Fig. ^[ but for GaAs spacers. Here, the 
values of 2n/(2k F ) are 40ML, 64.5ML and 95.2ML for E = 
20, 10 and 5mRy, respectively. 



deed, minority-spin states incident from Fe at Ep having 
the A 5 (C , 4 1 ,) and A 2 (C4„) symmetry (see Fig. [l]) cannot 
couple to the SC Ai(C2„) conduction band, but they 
can couple to decaying SC states that have the correct 
symmetry. In this way, if the thickness of the spacer is 
moderate, they can have an important contribution to 
the current through tunneling.c3 For larger thicknesses 
they become unimportant, and the asymptotic oscilla- 
tory behaviour appears. This situation of co-existence of 
tunneling current with "normal" current is much stronger 
in GaAs, because it has a smaller band gap than ZnSc, 
and thus the decay length of such evanescent states is 
much longer. 

In Figs. | and |, the majority-spin <?(ky) is demon- 
strated for 97ML-thick spacers of GaAs and ZnSe, for 
gate voltage shifts of Eq — 5mRy, lOmRy and 20mRy. 
The conductance resonances form rings around k|| = 0, 
up to hp; they are what one expects by rotating the 
graphs of Fig. || around the origin. It is remarkable that 
the majority-spin conductance is quite isotropic in all 
cases. In contrast, we find the minority-spin conduc- 
tance rings to reflect more the quadruplicate structure 
of the surface Brillouin zone, but seem actually (by in- 
spection) to obey one extra symmetry operation and to 
be octuple, as if the group were C^ v . This is observed in 
all cases, and is mostly evident in the case of ZnSe with 
Eq = 20mRy, where kp is largest; this is shown in Fig. |t]. 
We shall give the explanation of these observations to- 
gether with the analysis of similar data for the antiparal- 
lel magnetic configuration of the leads, at the end of the 



Fe / 97ML GaAs / Fe Majority 



1.0 
0.5 



CD 

CD 
O 
C 

o 

3 
"D 
C 

o 
O 



1.0 
0.5 




k p =0.050 
20mRy 



k F =0.031 
10mRy 



1.0 
0.5 



k p =0.021 
5mRy 



FIG. 5: Conductance (ky-resolved) of majority-spin electrons, 
in the case of a Fe/ 97ML GaAs /Fe junction. Top: E = E c + 
20mRy, k F = 0.050; Middle: E = E c + lOmRy, k F = 0.031; 
Bottom: E — iJc + 5mRy, k F = 0.021. The ky axes are along 
the f — M directions. 



next section. Evidently, the majority-spin conductance 
retains its dominance over its minority counterpart; the 
kii-integrated conductance is presented in Tabic |[ 



VII. ANTIPARALLEL MOMENT IN THE 
LEADS - MAGNETORESISTANCE 



From the analysis presented in Sections [V and [v| one 
should expect a strong reduction of the conductance if 
the magnetic moments of the leads have an antiparallel 
orientation. If the moment of, say, the second lead is 
reversed, then the majority and minority bands will be 
interchanged there. So, the incoming minority-spin elec- 
trons will be nearly blocked at the first interface, while 
the incoming majority-spin electrons will propagate up 
to the second interface but suffer almost total reflec- 
tion there, since they will encounter the states of the 
A2'(C4„)-type to which they do not couple well. Again 
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FIG. 6: Conductance (k||-resolved) of majority-spin electrons, 
in the case of a Fe/ 97ML ZnSe /Fe junction. Top: E = E c + 
20mRy, k F = 0.083; Middle: E = E c + lOmRy, k F = 0.056; 
Bottom: E = E c + 5mRy, k F = 0.038. The kn axes are along 
the r — M directions. 
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FIG. 7: Conductance (ky -resolved) of minority-spin electrons, 
in the case of a Fe/ 97ML ZnSe /Fe junction for E = _B C + 
20mRy, &f = 0.083; an octuple symmetry is evident. The kii 
axes are along the T — M directions. 
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FIG. 8: Fe/GaAs/Fe (top) and Fe/ZnSe/Fe (bottom) conduc- 
tance per spin channel along T — M (k x ), for the antiparallel 
magnetic configuration of the leads, and for several spacer 
thicknesses. The values are much lower than the majority- 
spin conductance and much higher than the minority-spin 
conductance in the parallel case. The fef-values are: 0.021, 
0.031, and 0.050 for E = 5, 10, and 20mRy, respectively. 



the situation is analogouS|-to the one encountered in the 
case of tunneling barriers .t3 

Indeed, in Fig. [s] we see the that the conductance in 
the antiparallel configuration is calculated to be orders 
of magnitude lower than the majority-spin conductance 
(and the total one) of the parallel configuration, but still 
orders of magnitude higher than the minority-spin con- 
ductance of the parallel configuration. The effect can 
be understood in terms of the reflectance and transmit- 
tance at the interfaces. If Tj- is the (high) single-interface 
transmision probability involving majority Fe states and 
T s \ is the (low) one involving minority Fe states, with 
case of parallel alignment the ma- 
jority electrons will have a total transmission probability 
from both interfaces of the order of T to t ~ (T^) 2 (neglect- 
ing resonance effects), the minority ones (T^) 2 , while the 
antiparallel-configuration electrons will have T s ;T S j for 
each spin channel. Evidently, (T s \) 2 > tTtJ- > (T^) 2 , 
q.e.d.. We observe, by the way, that this line of thought 
suggests that in the antiparallel configuration the con- 
ductance gn (per spin channel) is the geometrical aver- 
age of the conductances of the two spin channels in the 
parallel case: gn — y/gYtWll (to be valid but for backseat- 
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FIG. 9: Conductance (ky -resolved) of incoming majority elec- 
trons, in the case of a Fe/ 97ML GaAs /Fe junction with 
antiparallel magnetic orientation of the two Fe leads, for 
E = E c + 5mRy, fci? = 0.021; a quadruplicate symmetry 
is evident. The ky axes are along the F — M directions. 



tering effects). This is true for k|| in certain directions of 
the surface Brillouin zone, i.e. along k x and k y (the cu- 
bic axes), including of course k|| = 0. At such fc-points, 
the transmission through the first interface (Fe into SC) 
is the same as through the second (SC into Fe); how- 
ever, for other ky -points this is not true, so spin-up and 
spin-down electrons have different <?(kii) and only equal 
ky -integrated g as shall be explained in the end of the 
section. Our numerical results verify this. So, for an ar- 
bitary k||-point, the geometric average relation can hold 
at most for the order of magnitude. We note in pass- 
ing that, if we had a spacer material with C^ v interface 
symmetry, as e.g. MgO, the geometric average rule would 
not hold at all, because the minority Fe A2' (C^-state 
would be orthogonal to the spacer A^C^) conduction 
band; then the minority electrons would reach the sec- 
ond interface only through a complex band with expo- 
nentially damped probability and the assumptions of the 
two-reflectance-argument would not hold. 

For a large spacer thickness of 97MLs we see in Fig. ^ 
the ky -resolved conductance for the transmission from in- 
coming majority-spin to outgoing minority-spin channels. 
In Table |, we see the integrated (over the SBZ) conduc- 
tance for several gate voltage parameters Eq in the case 
of 97ML-thick spacers, together with the spin current po- 
larisation P — (gn — 9ix)/(ffTT + .9ll) an d the magnetore- 
sistance (MR) ratio defined as (g^ +.9ti)/ (.9TT + .9n) (* ne 
so-called "pessimistic definition"). Evidently the calcu- 
lated device acts as an almost ideal spin filter and switch 
with extremely high MR ratio. For lower energy shifts 
the spin filtering and MR ratio increase, because the al- 
lowed k|| close up to f and the states have more and 
more Ai(C2,j)-charactcr. Because of the A2' minority- 
spin state, however, the ideal 100% cannot be reached 
even in the limiting case; in contrast, it would be reached 
e.g. in thepase of an MgO spacer because it exhibits C^ v - 
symmetry.e3 

As promised at the end of the previous section, we 



now turn our attention to the explanation of the circu- 
larly symmetric form of g(k||) for the majority electrons 
in the parallel-alignment case, vs. the octuple symmetry 
seen for the minority electrons, and all this vs. the qua- 
druplicate symmetry in the antiparallel-alignmcnt case. 
As ky departs from f , the Fe and SC states do not 
belong exclusively to a single representation any more, 
but are rather admixtures of the various representations; 
but they still retain mostly the character they had at 
I\ In the language of localised orbitals, the majority- 
spin states are formed mostly by the circularly symmetric 
s+p z +d z 2 orbitals (plus small admixtures away from T); 
the minority-spin states consist of d xy from the A2' (C^ v )- 
band, p x + p y + d xz + d yz from the A5(C4i,)-band, and 
d x 2_ y 2 from the A2(C4„)-band; finally the SC conduction 
band states consist of s + p z + d xy from the Ai(C2i,)- 
band. Away from T, new orbitals start to contribute to 
each band, but in amounts negligible for our discussion, 
since we remain close to T. 

Firstly we concentrate on the coupling of the minority- 
spin states. At exactly T, the only combination that gives 
nonzero inner product is d xy orbitals of Fe with d xy - 
like states of the SC; the rest of the combinations are 
inner products of symmetric with antisymmetric wave- 
functions, resulting to zero. As kii departs from T, the 
p x , p y , d xz , d yz , and d x 2_ y 2 minority states of Fe atoms 
neighbouring a particular SC atom at the interface ob- 
tain slightly position-dependent phases as e lk n r ; then 
the wavefunctions formed by combining them obtain a 
small part symmetric around the SC atom, and this gives 
nonzero inner product with the SC s -t- p z -t- d xy . This 
overlap integral, in first approximation proportional to 
k\\ , is different for the various directions of kii , following 
the pattern of d xy . Clearly then the bonding and the 
conductance must have a quadruplicate symmetry in ky- 
space, as does d xy in real space. By inspection of the 
Fig. [7] we see an octuple symmetry. The explanation for 
the extra symmetry lies in the zincblende geometry and 
the directionality of the bonding. Indeed, as we enter the 
SC (e.g. ZnSe) from the one lead, we encounter Zn and 
then Se on the tetrahedral positions along the (x, y) di- 
agonal; but as we leave it, we encounter Se and then Zn 
on the tetrahedral positions along the (x, — y)-diagonal. 
Thus, the directionality of the SC d^-like states and con- 
sequently the bonding and transmission properties of the 
two interfaces are equivalent but rotated by 90° to each 
other, so the combined transmission obeys one extra sym- 
metry operation and is octuple. 

Secondly we focus on the coupling of the majority- 
spin states. There the situation is simpler: Fe has only 
s+Pz+d z 2 circularly symmetric orbitals which can couple 
only to the SC s +p z , but not to d xy . Thus no direction- 
ality is induced by the latter; even as we depart from T, 
the small difference in phase obtained by neighbouring 
Fe sites gives only an antisymmetric part to the com- 
bined wavefunction and this has still zero inner product 
with the SC s + p z + d xy . The result is that the bond- 
ing and transmission properties for majority are isotropic 
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TABLE I: Calculated current polarisation and magnetoresistance (MR) ratio in the case of 97ML-thick spacers of ZnSe and 
GaAs for several gate voltage shifts Ef — E c ; both are close to the ideal 100%. The spin-dependent conductance g integrated 
over the surface Brillouin zone is also shown for both cases of magnetic orientation of the leads: parallel (majority and minority) 
and antiparallel per spin (the same for the two spin channels). 



Material 


Ef — E c 






g(e jh) (per umt-ccll surface area) 




Polarization 


MR ratio 






Majority 




Minority 


Antiparallel /spin 








5mRy (68meV) 


1.6 x 10" 


7 


1.0 x 10" 12 


2.6 x 10" 1() 


99.999% 


99.678% 


GaAs 


lOmRy (136meV) 


7.1 x 10" 


7 


2.1 x 10 -11 


2.7 x 10 -9 


99.994% 


99.229% 




20mRy (272meV) 


1.9 x 10" 


6 


5.3 x 10" 11 


7.8 x 10" 9 


99.994% 


99.196% 




5mRy (68meV) 


1.7 x 10" 


7 


2.4 x 10" 12 


2.0 x 10" 9 


99.972% 


97.746% 


ZnSe 


lOmRy (136meV) 


8.2 x 10" 


7 


3.0 x 10" 10 


1.4 x 10" 8 


99.926% 


96.553% 




20mRy (272meV) 


2.8 x 10" 


6 


2.4 x 10 -9 


6.7 x 10~ 8 


99.823% 


95.128% 



arround T. 

Finally we look at the antiparallel magnetic configura- 
tion of the leads. There, one either enters with circularly 
symmetric transmission via majority and exits with the 
quadruplicate symmetry via minority with a quadrupli- 
cate net result, as seen in Fig. or, for the opposite spin, 
enters with quadruplicate symmetry via minority and ex- 
its with circularly symmetric transmission via majority, 
again with a quadruplicate net result. For the two last 
cases, by the way, the gf(kii ) are rotated to each other by 
90°, again due to the aforementioned direction difference 
in the bonding; thus only along k x and k v is <?(k||) the 
same for the two spin directions in the antiparallel case. 

The same symmetry of S , (ky) as hexe is seen in results 
for tunneling Fe/ZnSe/Fe junctions^! so once more we 
see the formal connection between spin injection and tun- 
neling. 



VIII. LIMITATIONS AND SUMMARY 

Before summarising, we shall briefly discuss the limita- 
tions of our approach and the relevance to realistic exper- 
imental situations. Two main points must be addressed 
here: (i) the influence of diffuse, k|| -violating scattering 
and (ii) the possible effect of a Schottky barrier. As for 
point (i), it is true that the formation of terraces or steps 
and interdiffusion lead to diffuse scattering. To what ex- 
tent this reduces the control over the conductance must 
be examined seperately in each case and is a huge but 
challenging task. In the case of the Fe/GaAs interface it 
is known that in growing of Fe on GaAs the As atoms act 
like surfactants forming always an As monolayer on Fe. 
This is of course an indication that the interface struc- 
ture is not perfect. But progress is being done and one 
can reasonably hope that the quality of the interfaces will 
increase a lot in the future. 

About point (ii), Schottky barriers are known to ex- 
tend over mesoscopic lengths, especially when the doping 
is low. However, techniques to use quantum well struc- 
tures have resulted in lowering the conducion band un- 



der the Fermi level without direct impurity doping; such 
a situation would be modeled by our "gate voltage" pa- 
rameter Eq in addition to a real gate voltage. Then the 
Schottky barrier would be much shorter, in fact being 
determined by the Fermi level pinning due to the metal- 
induced gap states. On the other hand, in the single- 
interface calculations for spin injection by Wunnicke et 
al.c3 the effect of a Schottky barrier has been studied by 
emulating it with a long region near the interface where 
the SC potentials were kept to their physical unshifted 
positions, and the electrons had to really tunnel into the 
conduction band. The result was quite encouraging, giv- 
ing still an extremely high current spin polarisation. 

To summarise, we have performed ab initio cal- 
culations of the spin-dependent transport through 
Fe/GaAs/Fe and Fe/ZnSe/Fe (001) junctions, with a 
gate voltage parameter acting on the semiconductor so 
that the Fermi level lies slightly in the conduction band. 
The electron transport was supposed to be completely 
ballistic, assuming a perfect interface structure and two- 
dimensional periodicity perpendicular to the direction of 
growth. Under these assumptions we have shown that 
such systems can exhibit an extremely high degree of 
current spin polarisation and also a magnetoresistance 
ratio approaching the ideal 100%. We have been able 
to trace down these nice properties to the difference in 
the bulk band structure for the two spin directions of 
Fe, and also to the difference in the bonding of majority- 
and minority-spin states at the interface with the semi- 
conductor. In the same terms we have explained the high 
magnetoresistance values. We have also examined inter- 
esting interference effects that show up in such a junction 
due to the presence of two, rather than one, interfaces, 
and discussed the question whether these effects can in- 
vert the detected current polarisation. 

We have seen that the understanding of these systems 
stands in close connection with the understanding of bal- 
listic magnetic tunnel junctions, if one formally replaces 
the band structure near the center of the conduction band 
of the semiconductor with the complex band structure 
in the gap region. In both cases, it is important that 
very few states perform the conduction, namely the ones 
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near the center of the surface Brillouin zone; to know the 
properties of these states means to have control over the 
conductance. 

We have concluded that the control over the desired 
properties of such systems is best when one deals with 
ballistic transport. Diffuse scattering, particularly at the 
interface, would intermix the various conducting channels 
and cause the injection efficiency and magnetoresistance 
to drop; on the other hand, clean and abrupt interfaces 
preserve k|| and act as spin-selective transmitters and 



detectors. 
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43 Similar interference resonances in g(k\\) haus been ob- 
served in calculations of tunneling junctions) 4 -!! for trans- 
mission through those evanescent states that also have a 
nonzero real k z . Of course there the transmission probabil- 
ity is always many ordersrefrmagnitude smaller than one. 

44 Recent theoretical articlesEj'ES based on model calculations 
have suggested that a modulation of the current polarisa- 
tion due to interference by variation of &f is possible. Our 
arguments deny this possibility for the systems considered 
here, but surely not for all systems. Note that Refs. ^ 
and |3^ refer to systems of lower dimensionality (2D or 
ID), where the averaging over kn has a lesser effect (or is 
nonexistent, in ID). Also, the aforementioned modulation 
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would be perhaps possible also in 3D if the minority-spin 
reflectance of the interface were lower. 
The situation is analogous to the one-dimensional case of 
a free electron of wavenumber fc n = \[E encountering an 
infinitely long step barrier of height Vb where the wavenum- 
ber is fcout = VE — Vq. The transmition probability there 
is T s i = 4fc ou tfcm/(fcout + fc in ) 2 going to zero as y/E — Vq 
for small fc ou t- The point is that the group velocity in 
the barrier region goes to zero, although the transmission 
amplitude t remains finite, so current conservation forces 
T s i = \t\ 2 k out / ki n also to go to zero. The same point is trpe 
in our case when we approach the conduction band edgepj 
and thus the analogy is valid. 

This applies also for the k\\ -dependence of g: one sees the 
multi-resonant form of Fig. |^ if one substitutes the value 
for T s i(fc|| = 0) taken from the single interface calcula- 
tion into eq. (113), and takes for the other ky's T s i(k||) = 

T a i (0)k z /kp = r s i(0)«/fe|, — k?,/kjp; in such an approxima- 



tion, of course, any anisotropy of T S i(kii) is lost. 
This small imaginary part e, in principle infinitesimal, can 
also be viewed as imposing an attenuation exp(— y/eD) on 
the waves crossing a spacer of thickness D. Its effect can 
be studied by replacing to. — > £2exp( — \/tD) and r-i — * 
V2 exp(—2y/fiD) in eq. (|10[). It turns out that the resonance 
maxima are then reduced for larger thicknesses, and that 
the reduction is stronger for smaller T^. Thus the minority- 
spin conduction maxima are much more affected. 
Actually, the maxima in the DOS do not always have to co- 
incide with those in the transmission; a relative phase shift 
can occur, as pointed out in Ref. However, such a shift 
is not expected in the present case if we follow the analysis 
of Ref. fid, since there are no zeros of the transmission. 
The situation bears again an analogy with the magnetic 
tunnel junctions; in the limit of large spacer thickness, the 
currepti polarisation and MR ratio dortjot-ireach 100% for 
ZnSejB but they do reach it for MgO.BB 



